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Abstract 

Some of the classical orthogonal polynomials such as Hermite, Laguerre, Char- 
lier, etc. have been shown to be the generating polynomials for certain combinatorial 
objects. These combinatorial interpretations are used to prove new identities and 
generating functions involving these polynomials. In this paper we apply Foata's 
approach to generating functions for the Hermite polynomials to obtain a triple 
lacunary generating function. We define renormalized Hermite polynomials hn{u) 
by 

oo „ 



EK{u)— = e" 



n=0 



and give a combinatorial proof of the following generating function: 

2^ ^n[u)^^ - 2^ 23"(3n)!(l - 6u;z)3" (2n)! ' 

71=0 n=0 



where w = {I - VI - l'2.uz)/Qz = uC{Zuz) and C{x) = (1 - VI -4x)/(2x) is 
the Catalan generating function. We also give an umbral proof of this generating 
function. 



1. Introduction 

The Hermite polynomials Hn{u) may be defined by the exponential generating function 

Y.Hn{u)'- = e^-^\ (1) 



n=0 
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In this paper we will give a combinatorial proof of an identity for Hermite polynomials. For 
our combinatorial interpretation, it is more convenient to take a different normalization 
of the Hermite polynomials, which makes all the coefficients positive. Therefore, we work 

i'^ f —iu\ 



with the polynomials hn{u) = ^^Hn y—j=-j , where i = y— which have the generating 
function 

n=0 

All of our formulas for hn{u) are easily converted into formulas for Hn{u). 

Foata gave a combinatorial proof of Doetsch's identity [2J giving a generating 
function for h2n{u): 

EM«)^ = (l-2.)-/^exp ^^ . (2) 

n=0 ■ ^ ^ 

We will prove the following generating function for h^n{u): 

2^/i3n(M)^ = /i_gy^^-Z^23-(3n)!(l-6wz)3"(2n)!' 

n=0 * n=0 \ / \ / \ / 

in which w = {1 — a/1 — Vluz)lQz = uC{3uz), where C{x) = (1 — — 4x)/2x is the 
Catalan number generating function. Note that the formula can be written in terms of 
hypergeometric series: 

°° " Jw~u){3u~w)/6 / g 54^^ 



^ " n! a/1 — Qwz \6'6' '{l — QvjzY 



where 2^o(a, z) = y^(a)n {b)n—;, and 

n>0 

{a)n = a{a + 1) . . . (a + n — 1) 

is the rising factorial. 

We prove formula Q by two methods — umbral and combinatorial. In section |21 we 
define an umbra and study some of its properties and then give the umbral proof. An 
umbral proof of a generating function for h2m+n{u) is given in p. 

In section |3] we prove Q combinatorially by showing that both sides enumerate the 
same weighted objects. We first describe the combinatorial interpretation of the Hermite 
polynomials and then give the details of the weighted objects counted by both sides of 
the formula. 

By using these methods, it would be possible to give umbral and combinatorial proofs 
of a more general generating function for h3m+2n+k{u) . 
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2. The Umbral Proof 



Rota and Taylor in [7j laid a rigorous foundation for the classical umbral calculus. They 
consider a vector space of polynomials in several variables or "umbrae" and define the 
linear functional eval on it. A sequence (a„) is represented by an umbra A if evaA{A"') = a„ 
for all n. In practice, the word eval is usually dropped and we simply write = an with 
the understanding that the functional has been applied. When we write f{A) = g{A), we 
mean eval(/(y4)) = eval{g{A)). We consider formal power series f{t) with coefficients in 
a ring of formal power series R = Q[[x, y, z,- ■ ■]]. 

oo 

Definition 2.1 A formal power series f{t) = fnt^ 'is admissible if for every monomial 

n=0 

x'^yiz^ ... in R, the coefficient of x'^y^ z^ ■ ■ ■ in fn is nonzero for only finitely many values 
ofn. 

So, for example, f{t) = e^* is an admissible formal power series, while f{t) = e* is not. 
Some computations similar to those of this section can be found in section 4 of Gessel 



2.1. The umbra M and its properties 

We now define an umbra M whose relation to the Hermite polynomials will be described 
in the next section. In this section, we study some properties of this umbra that we will 
need. We define the umbra M by e*^^ = so that 

M" = <j 2''k\ (4) 
I 0, if n is odd. 

The following two formulas hold for M. 
Lemma 2.2 (i) For any admissible formal power series f, e^^'f{M) = e''^^f{M + z). 



Proof. (i) We first observe that it is sufficient to prove that the formula holds for 
/(t) = e*^': If the formula is true for f{t) = e*^, then comparing coefficients of x^/n! on 
both sides we get that the formula holds for f{t) = t". But this implies by linearity that 
the formula is true for all admissible formal power series /. To prove the formula for 

f{t) = e*^, we have 

Mz Mx _ M{z+x) _ z'^/2+zx+x'^ /2 _ z'^/2 zx Mx _ z'^/2 (M+z)x 

(n) We have 

M-^z ^,.2n^" ^{2n)\z- 1 



e 



^ n\ ^ 2"ra! n! J\ - 2z 

n=0 n=0 



□ 
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Corollary 2.3 For any admissible formal power series f, 



e^'^fiM) = -=!== f( 



Proof. As in the lemma, it is sufficient to prove the formula for f{t) = e*^. Here 
gM = z+Mx ^ apply Lemma f2. 21 ii) directly, we cannot eliminate the linear 

term in M. So we introduce a parameter a and rewrite z+Mx ^Ma^M z+M(x-a)_ 
will choose a value for a later. Now applying Lemma 12.21 (i) we get 

Ma M'^z+M{x-a) _ a'^ /2 {M+a)'^z+{M+a){x-a) _ M^z+M{x-{l-2z)a)+za'^+xa-a^ /2 

Now we choose the value of a to eliminate the term in M on the right; i.e., we solve 
X — (1 — 2z)a = and get a = x/{l — 2z). Substituting this value of a in the above 
expression and simplifying, we obtain that z+Mx _ z^x /2(i-2z) ^ gy Lemma [ 

this is equal to g^V2(i-2g)^ applying the definition of M directly gives 

a/I - 2z 

1 J M 



— ] ^ pMx/y/l^ ^ ^x^/2{l-2z) 



VI - 2z VVi - 2^7 VI - 22 VT^ 



Applying the corollary to f{M) = e^^^^ gives the formula 



□ 



exp 



2^)3/2 



vr^2i 

which we will need in the next section. 



(5) 



2.2. Proof of the formula 

To prove formula Q, we first express the Hermite polynomial /i„(m) in terms of the umbra 
M. We have 

OO ^ 

yh{u^^ = eUZ+z^/2^^{u+M)z_ 

Comparing the coefficients of on both sides, we get hn{u) = {u + M)^. Using this, 

we get 

OO OO ^ 

y h,Ju)^ = y(u + Mf^^— = e("+^^)'^ = e'^n^+zwiu^+zM^u+Nfi)z_ 

n=0 n=0 
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We follow the same procedure as in the proof of Corollary 12.31 We introduce a parameter 
a and rewrite the last expression as e^"e"^^+*^^^"^^~")+^^^^"^"^^^^^. Now applying Lemma 
O (i) we get 

_ gZQ3+3n^a2-a2/2+3«^za+M^z+Af{3za^+(6u2-l)a+3n2^)+M2(3zQ!+3«z)+A/^z ^g-j 

In order to apply formula (0), we need to eliminate the linear term in M . So we choose 
a value of a that makes the coefficient of M equal to zero. By solving a quadratic and 
taking the solution with a power series expansion, we get a = (1 — \/\ — 12uz)/Qz — u. 
We know that the Catalan generating function C{x) is given by (1 - VI - 4x)/2a;. In 
terms of this generating function, a = uC{3uz) — u = w — u where w = uC{3uz). Using 
this value of a to simplify expression we get that 

°° 

hn(u)— = expiw^z ~ (w- ufl2 + BM^wz + M^z). 

n=0 

We know that C{x) satisfies the equation C{x) = 1 +x(C(x))^. Substituting x = 3uz 
and using the fact that C{3uz) = w/u, we obtain w = u + 3w'^z] i.e., w — u = 3w^z. This 
gives us 

f2,,,_ 1 » ,r3 



^{w-u){3u-w) /6^3APwz+M'>z 



y^hniu)— 

n=0 

Applying formula (0) with f(t) = t^z, we get that 

- n g{^-«)(3«-»)/6 / ^3^ X 



Then writing the second exponential function as a series and using (^, we obtain the final 
result: 

^{w-u){3u-w)/e (6n)! 



n\ VI - Qwz ^23"(3n)!(l-6wz)3'^(2n)!' 

n=0 n=0 \ / \ y \ y 

Now we turn to the combinatorial method of proof. We begin with a combinatorial 
interpretation of the Hermite polynomials that will be used in the combinatorial proofs. 



3. The Combinatorial Proof 

We assume that the reader is familiar with enumerative applications of exponential gen- 
erating functions, as described, for example, in Chapter 5] and The product 
formula and the exponential formula for exponential generating functions discussed in 
these references play an important role in the combinatorial proofs. 
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3.1. Combinatorial interpretation of Hermite polynomials 

The exponential generating function uz + z^ /2 counts sets with one or two elements, where 
a one-element set is weighted u. Then by the exponential formula, the coefficient of /n\ 
in e"^+^ which is hn{u), is the generating polynomial for partitions of an n-element 
set into blocks of size one or two, where each block of size one is weighted u. (If we 
used Hn{u) as Foata did, instead of hn{u), we would need to attach a weight of —2 to 
each two-element block and a weight of 2u to each one element block.) It is convenient 
to represent these partitions as graphs in which the vertices in a two-element block are 
joined by an edge. We call these graphs, in which every vertex has degree at most one, 
matchings . 

Thus the Hermite polynomial hn{u) can be viewed as the generating polynomial for 
the number of vertices of degree zero over the set of all matchings on n vertices, where 
each vertex of degree zero is assigned the weight u. With this combinatorial interpretation 
we will prove the following formula: 

V h-, (U)— = p(^-n)(3n-^)/6 ^ ^ (^^ 

^ '"^ ^n! J\ - Qwz ^ 23«(3n)! (1 - ^wzf^ (2n)!' ^ ' 

We will describe the graphs enumerated by the left side of formula (|7j). Then we will 
describe the same graphs in terms of their connected components and use the product 
formula for exponential generating functions to complete the proof. 

3.2. Graphs counted by the left side. 

In order to give a combinatorial interpretation to Yl'^=o^3niu)x"'/n\, we must interpret 
hsn{u), which counts matchings of a 3n-element set, as counting labeled objects with n 
labels. To accomplish this, we take n labels and attach to each one three vertices marked 
a, b, and c. Figure ^ shows a labeled vertex connected to three marked vertices. Then 



c 

Q 




a 



b 



Figure 1: A label and three marked vertices 

h^niu) counts graphs constructed by taking n components like that in Figure [H with 
labels 1, 2, . . . , n and adding to them a matching of the 3n marked vertices, where each 
unmatched marked vertex is weighted u. Figure |21 shows such a graph. 

Let G be the set of all graphs enumerated by the left side of formula ((Zj). For a graph 
in G, each vertex has degree one, two, or three. The trivalent vertices are labeled with 
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the integers 1 to n, where n is the total number of trivalent vertices, and are weighted by 
the exponential generating function variable z. The bivalent and monovalent vertices are 
marked with a, b, and c; the monovalent vertices have weight u and the bivalent vertices 
have weight 1. 

We first make a preliminary simplification of the graphs in G: we eliminate all the 
bivalent marked vertices, moving their "marks" to the adjacent trivalent vertices. 

More precisely, we think of each edge joining two trivalent vertices as consisting of 
two "half-edges" , each of which has a mark. Although we retain the monovalent vertices, 
we move their marks to the half-edge of the adjacent trivalent vertices. Note that some 
trivalent vertices now have loops. Figure which shows the transformed version of the 
graph in Figure |21 should make this simplification clear. 

We will express the generating function for graphs in G as a product of three factors 
corresponding to connected components with no cycles, with exactly one cycle, and with 
at least two cycles. Note that a loop is a cycle. 

3.3. Counting w-trees 

The graph in Figure 121 has three connected components. One component is a tree. The 
component with one cycle may be viewed as a cycle of trees. The third component may 
also be decomposed into trees and cycles. Thus the first step of our enumeration is to 
count the rooted trees from which we will construct our graphs, which we call w-trees. 

Definition 3.1 A w-tree is a rooted tree with labeled trivalent internal vertices and un- 
labeled leaves in which the three half- edges incident with every internal vertex are marked 
with the letters a, b and c, and one of the half- edges incident with the root has no matching 
half-edge. 

Figure m shows a w-tree. A special case of a w-tree is the one with no internal vertices 
and no marks: it is simply a leaf with a half-edge attached to it. 

We denote by w the exponential generating function in z for w-trees where each 
internal vertex is weighted with z and each leaf is weighted with u. Let Wn be the number 
of w-trees with n internal vertices. There are 2" ways to draw such a tree (with the root 
at the top) by choosing an ordering of the two children of each internal vertex, so there 
are 2"'Wn such drawings. But these drawings can be counted in another way. If we remove 
the labels and marks from such a drawing, we have a binary tree counted by the Catalan 
number C„ = (^^) . The marks can be added in 6" ways and the labels in n\ ways. Thus 
= 6"'n\ Cn, so Wn = 3"n! C„. Each such tree has n -\- 1 leaves, so the exponential 
generating function for weighted w-trees is w = Yl^=o S"'^' CnU'^'^^z^ /n\ = uC{3uz). This 
relation to the Catalan generating function gives the equation w = u + Sw'^z as shown in 
section 12.21 

3.4. Acyclic components. 

The connected components of graphs in G with no cycles are trees. Such a tree appears 
in the lower left corner of Figure El 
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Figure 4: w-tiee 



To count these trees we use a known method (see, e.g., [3]) for relating labeled unrooted 
trees of various types to rooted trees: we subtract the edge-rooted versions from the 
vertex-rooted versions. 

First we count trees rooted at a trivalent vertex. In such a tree, the three half-edges 
at the root are joined to w-tiees. So we can construct such a rooted tree by taking an 
ordered triple of w-tiees, with exponential generating function w^, and a new root vertex 
with its three half-edges marked a, b, and c, and attaching the first w-tiee to half-edge a, 
the second to half-edge b, and the third to half-edge c. Thus the exponential generating 
function for these vertex-rooted trees is w^z. 

Next we count versions of these trees rooted at an edge joining two trivalent vertices. 
With the help of Figure El we see that the exponential generating function for such edge- 
rooted trees is (3w^z)^/2 = 9w'^z'^/2. 



w-tree 



w-tree 



edge - root 




3 choices 
for this mark 




w-tree 



3 choices 
for this mark 



w-tree 



Figure 5: Edge-rooted tree 



An unrooted tree with n trivalent vertices is counted n times in the generating function 
for vertex-rooted trees and n — 1 times in the generating function for edge-rooted trees. 
Thus the difference w^z — 9w'^z'^/2 counts every unrooted tree once. A straightforward 
computation shows that w^z — 9w'^zy2 can be expressed most simply as {w — u){3u — w)/6. 
Thus the exponential generating function for graphs whose connected components are 
trees is 
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As pointed out by the referee, other derivations of this generating function can be 
obtained by using the combinatorial interpretation of derivatives. Here is one such deriva- 
tion: If T is the generating function for these trees, then 



du 



w — u 



1 - VI - 12u2 
6z 



— u 



(9) 



since if we remove one of the monovalent vertices (weighted u) from a tree counted by 
T, leaving an unmatched half-edge, we obtain a w-tree with at least one trivalent vertex. 
Integrating with respect to u, and using the fact that T has constant term 0, we get 

(1 - Vluzfl^ 



T 



u 



u 



108^2 + 2 ' 



and this is easily checked to be equal to {w — u){3u — w)/6. 
We also note that the coefficients of T are given explicitly by 



n=l 



(2n)! 
(n + 2)! 



u 



n+2' 



Now let us look at the connected components with one or more cycles. The two 
components other than the tree in Figure El are of this form. First we reduce these 
components by shrinking all the w-trees present in them. Figure IHl illustrates this process. 
As we see in the picture, the process of shrinking w-trees leaves behind vertices weighted 




Figure 6: Components with cycles with w-trees eliminated 



with z which are now bivalent. We further reduce the components by eliminating these 
newly created bivalent vertices. Figure [7| shows what we get after this further reduction. 
Thus we get either circles with no vertices or components (with multiple edges and loops) 
in which each vertex is trivalent. We will consider the two cases separately. 
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Figure 7: Graph in Figure IHl with bivalent vertices ehminated 



3.5. Components with one cycle. 

If the reduction of a component of a graph in G results in a circle, then the original 
component has exactly one cycle and each trivalent vertex on the cycle has a w-tree 
attached at its third half-edge. The exponential generating function for directed cycles is 
z'^ In. To count directed cycles with w-trees attached, we replace z by wz, and to 
add three marks to the half-edges of each vertex in the cycle, we multiply {wz)'"' by 6"". 
Finally, to undirect the cycles, we divide by 2. (Note that marking the half-edges for the 
cases n = 1 and 2 destroys the symmetry that would prevent us from dividing by 2 before 
the marks are added.) 

Thus the exponential generating function for cycles of w-trees is 



and so the exponential generating function for graphs whose connected components are 
cycles of w-trees is 



3.6. Components with more than one cycle. 

The components of graphs in G with more than one cycle reduce to connected graphs with 
only trivalent vertices. The connected graph on the left in Figure [7| is of this type. Rather 
than counting connected graphs of this type, as we did in the previous two sections, and 
then exponentiating, we count graphs, not necessarily connected, whose connected com- 
ponents are of this type. These graphs with m labeled vertices are precisely the graphs 
counted by h^miu) with no monovalent vertices. To see this, recall that in section EISl 




n=l 



log (1 — 6ui2:) 



1 



(10) 



— 6wz 
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we described how each labeled vertex has three marked vertices attached to it and after 
constructing the matching on the marked vertices, the bivalent (or matched) marked ver- 
tices are eliminated as a preliminary simplification. If there are no unmatched marked 
vertices, the preliminary simplification gives a graph with only trivalent labeled vertices 
which have marks on their half-edges. Thus these graphs are counted by /t3m(0), the num- 
ber of complete matchings of 3m vertices. This number is for m odd and {6n)\/2^^{3n)\ 
for m = 2n. 

We can recover the original graph from the reduced graph by introducing an ordered 
sequence of trivalent vertices on each edge, giving each new vertex marked half-edges 
and attaching a w-tree to one of the half-edges. We can see this for the component with 
trivalent vertices in Figure[7|by tracing it back to its original component in FigureEl Thus 
we get a factor YlT=o (6^2;)'^ = 1/(1 — 6wz) for each of the 3n edges in the reduced graph, 
where the 6 is the number of ways to specify the marks at a trivalent vertex. Hence the 
exponential generating function for graphs whose components have more than one cycle 
is 

(Qn)] 1 



(11) 



„=o 2^"(3^)' (1 - 6wz)3" {2n)\' 

By the product formula for exponential generating functions, we multiply (jH)), ()1()|1 
and (fTTjl to get 



e 



(w—u){3u—w)/6 



E 



(6n)! z 



2n 



VI - Qwz ^ 23"(3n)!(l - Qwzf-^ {2n)\'' 



as the exponential generating function for all graphs in G, and this is the right side of (|7j). 
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